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However the lower limit for wave length of sound waves that 
can still be properly represented increases considerably : 



10 2 



(49) 



All numerical hydrodynamical codes that employ artificial vis- 
cosity in some form can therefore not model winds in which 
sound waves play a role. This essentially leaves implicit or high 
resolution shock capturing (HRSC) numerical schemes as the 
only option, since these do not need artificial viscosity to en- 
sure the numerical stability of shocks (see e.g. van Leer, 1979 ; 
Korevaar & van Leer, 1988). 

For HRSC hydrodynamical codes there are still the prob- 
lems of sampling in time and sampling in space. Although the 
time step is not limited by the CFL condition, it is limited be- 
cause the sound wave is only well sampled for n t points per 
period. 



At < 



sw 



(50) 



Here _P SW is the minimum period of the sound waves to be mod- 
elled. Because the period of the sound waves can be quite short 
compared to the period of the shock waves, it would take many 
time steps to calculate wind models. 

The spatial sampling problem is much more serious. The 
effects of discrete spatial sampling of the flow are sometimes 
also referred to as numerical viscosity. All such effects are an 
unavoidable consequence of the discretization of the fluid equa- 
tions as opposed to artificial viscosity which is added as an ex- 
plicit term in these equations. The finite spacing of grid points 
causes a sampling of the sound waves that is less accurate as 
the wavelength of the sound wave decreases. If the sound wave 
is considered to be sufficiently well sampled for a minimum of 
n x points per wavelength the maximum grid spacing for a given 
period of the sound waves that are to be modelled is : 



Ax < 



_P SW min (« sw ) 



where Ax is the spacing of the numerical grid and min (v sw ) is 
the minimum of the sound wave propagation speed in the rest 
frame of the star. Since the wind is likely to have at least one 
transition from subsonic to supersonic outflow at any time, this 
means that near that point the wavelength in the Eulerian frame 
decreases to 0. This implies that the required grid spacing would 
also decrease to zero. Thus the required resolution is always 
higher than can be achieved with a finite number of grid points. 
Therefore even an HRSC type adaptive grid method can not 
properly conserve wave action. 

All schemes that solve the fluid equations on a discrete mesh 
in space and/or time suffer from this effect of numerical vis- 
cosity. None of these schemes can therefore properly conserve 
wave action. The only solution is to include an equation like 
(34) into the scheme of equations to be discretized and solved 
simultaneously with the equations of conservation of mass, mo- 
mentum, and energy. 



6. Conclusions 

In this paper several points are made. The first point is that in 
stellar winds the interaction of sound waves with strong shocks 
generated by pulsation does not remove wave action density of 
the sound waves. It is shown that the sound wave action density 
increases upon the passage of a shock. This is relevant to the 
models of the winds of Mira variable stars and other pulsating 
stars that lose mass through a stellar wind. 

The effective damping length of the sound waves in the rest 
frame of the star can increase considerably because of the pe- 
riodic injections of wave action density through the passage of 
shocks. The combination of pulsation driven shocks and sound 
waves could therefore justify the large values of the dimension- 
less damping length used by Pijpers and Habing (1989). How- 
ever, the functional form of the wave action density differs from 
the one used by these authors and it is necessary to do similar 
work for the combination of sound waves and pulsation before 
the effectivity of this model can be discussed. 

The time dependent hydrodynamical simulations of the 
winds of Miras carried out to date do not and cannot properly 
take into account the dynamical effect of the presence of sound 
waves. A proper treatment requires the solution of the transport 
equation for wave action density jointly with the equations for 
mass, momentum and energy. If sound waves are generated in 
the surface convective layers of cool giants in much the same 
way as they are in the convective layer of the sun, their dynam- 
ical effect on the winds of cool giants should not be ignored. 
This means that the paradigm for winds of AGB stars, exem- 
plified by the work of Wood (1979), Bowen (1988), and Bowen 
& Willson (1991), needs to be adjusted. 
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this shock has a constant amplitude over a certain region in the 
wind not only the sign but the also the value of the logarithmic 
gradient of the sound wave amplitude that appears in equation 
(45) remains unaffected by the passage of the shock. 

From the analysis of the previous section it is clear that 
the wave action density and therefore the wave amplitude can 
increase quite rapidly with radius if (A) > (A*). It is clear 
from energy conservation arguments that the amplitude must 
decrease to zero towards infinity. Physically the damping length 
must become shorter as the wave amplitude increases so that at 
some point the effective damping length A e ft > 0. In princi- 
ple the force from the sound wave pressure can therefore point 
inward on average in a part of the wind but it cannot do so 
throughout the wind. This means that the presence of shocks in 
the winds of Mira variables does not preclude that sound waves 
enhance the mass loss and/or the outflow velocity of the wind. 

The relative importance of the sound waves and the pulsa- 
tionally generated shocks in generating the mass loss is hard to 
assess on the basis of this work alone. In one of the models of 
Bowen (1992) the mechanical power delivered at the piston for 
a 5 Mq model for driving at the fundamental mode period is 
~ 10 Lq which amounts to a fraction of the radiative luminos- 
ity of ~ 2.4 10~ 4 . Pijpers and Habing (1989) show some results 
for a sound wave driven wind for a star with the same hydro- 
static density scale height as the model of Bowen (Pijpers & 
Habing model 7). For sound wave driven winds this is the only 
relevant parameter, so the wind of each of these models can be 
compared. Bowen (1992) reports a quite low mass loss rate of 
< 10~ 8 Mq /yr whereas Pijpers and Habing (1989) show mass 
loss rates of up to 10~ 6 M /yr for a mechanical luminosity in 
sound waves of ~ 3 10~ 5 times the radiative luminosity. From 
this it would seem that the sound wave pressure alone is more 
efficient at generating mass loss than pulsational shocks alone 
for at least some stars. However this impression may be false 
for two reasons. 

Firstly : because of the non-linear interaction between sound 
waves and periodic strong shocks described in this paper it is 
not at all clear what the situation is if both are present. Sec- 
ondly : the velocity amplitude of the sound waves at the base 
of the wind is ~ 10% of the sound speed for these models. 
This might appear to put the sound wave pressure in the regime 
where the linearized approach of Pijpers & Habing is still valid. 
However, in the wind the wave action density increases con- 
siderably, and therefore so does the amplitude, before decreas- 
ing to zero again. The amplitudes become so large that non- 
linear effects are important. As pointed out by Koninx & Pij- 
pers (1992) this means that the results of the calculations done 
entirely in the linearized approach may not be reliable. 

From figure 2 of Pijpers & Koninx (1992) it is seen that at 
larger amplitudes the sound waves generate a larger force than 
would be expected from an extrapolation of the linearized for- 
mulation. If this effect is primarily of importance in the region 
of the wind beyond the critical point this has the consequence 
that for the same final velocity in the wind, the mass loss rate 
is lower than predicted in the linearized models. If this effect 
primarily plays a role below the critical point, the mass loss 



is effectively increased. In the absence of actual modelling of 
non-linear sound wave driven winds, with and without periodic 
shocks generated by pulsation, it is premature to claim that ei- 
ther of these mechanisms dominate in the determination of the 
mass loss of giant and supergiant stars. 



5. Implications for numerical modelling 

The above results are important not only for the models of 
sound wave driven winds. The results of the numerical mod- 
elling of pulsation driven mass loss of Bowen (1988) and Wood 
(1979) show that pulsation can drive a mass loss only in com- 
bination with the radiation pressure on dust formed in the stel- 
lar wind. However the effect of sound waves is not taken into 
account in their models. It is possible that the combination of 
propagating shocks generated by the pulsation and sound waves 
generated by the surface convection zone can drive the mass 
loss of red giants even without the driving of radiative forces 
on dust. 

It is important to realize at this point that sound waves are 
suppressed in all numerical hydrodynamical codes for reasons 
of numerical stability and also because of the finite spacing of 
grid points. All explicit hydrodynamical codes use artificial vis- 
cosity to stabilize shocks. If it is not used numerical inaccu- 
racies will shift the shock back and forth over distances small 
compared to the grid spacing which causes the shock to quickly 
disperse. This is exactly the same effect as described in this 
paper but instead of actual sound waves perturbing the shock 
the sound is produced artificially by 'numerical jiggling' of the 
shock front. A high artificial viscosity is introduced to damp 
these oscillations and thereby prevent the dispersal of shocks. 
The effect of this high viscosity is not only to damp the numeri- 
cal production of sound but it also suppresses real sound waves. 
The damping length A is artificially shortened and the dynam- 
ical effects of sound waves are destroyed. The usual recipe for 
artificial viscosity (Richtmyer & Morton, 1967) uses a coeffi- 
cient viscosity rj of the form : 



T) = l 2 p \V ■ v| 



(47) 



Here / is a numerical smoothing length. The Reynolds' number 
of sound waves with this artificial viscosity is : 



Re = 



pASv 

n 



A8i 



l 2 \V 



(48) 



The limiting wave length of sound waves below which the vis- 
cosity dominates the dynamics can be estimated roughly by set- 
ting Re = 1 . Away from shocks the local velocity divergence is 
probably of the order of the velocity amplitude 6v of the sound 
waves. Using this the minimum wavelength that can be repre- 
sented is of the order of the numerical smoothing length /. At the 
shocks the divergence of the velocity can be orders of magni- 
tude larger than the velocity amplitude of the sound waves. This 
is precisely the reason why artificial viscosity is introduced. 
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no work has been done on the damping of acoustic waves in 

the atmospheres and winds of cool giant stars. There is some 

work available on the excitation and damping of sound waves 

in homogeneous radiative media (cf. Hearn, 1972 ; Mihalas & 

Weibel Mihalas, 1983). Using the results from that work yields 

a rough estimate of the intrinsic damping length of between 10 

and 1000 times the intrinsic wave length of the wave. Using 

these results, waves with an intrinsic wavelength in the range 

of 0.1 to 0.001 times the static density scale height are mosj. 

likely to have an effective damping length approaching infin 7 ^ 

ity. This corresponds with intrinsic wave lengths A and period^ 

r of the acoustic waves of : on 

i — i 

A ~ 4 10 6 m - 4 10 9 m oo 

r ~ 15 min —Id on 

o 

Note that all of these estimates are quite uncertain since thef?: 
rely on analyses of wave damping that are not really appropriate, 
and simulations of shocks in the winds of cool giants with afi 
incomplete treatment of the radiative transfer. $ 
The original models of Pijpers & Habing (1989) for the 
winds of AGB stars match the observed properties remarkably 
well, but for quite large values of the damping length A. Al- 
though this has been considered as a weakness of those models, 
it may be that such large values arise precisely because of the 
interaction between the pulsation driven shocks and the sound 
waves. 



4.2. The gradient of the wave pressure tensor 

The presence of sound waves in a stellar wind is felt through 
the gradient of the pressure tensor and not through the pres- 
sure tensor itself. One should be careful when calculating the 
average value of this gradient over a pulsation period. Gradi- 
ent operators and the integral of the averaging operation do not 
necessarily commute because of the presence in the wind of 
the discontinuities from the pulsation. Koninx & Pijpers (1992) 
give an expression for the gradient of the non-linear wave pres- 
sure tensor for sound waves with a purely radial wave vector : 



(V • P) r 



-c]pl(l, 1 + l) 
1 du 1 



2 1 

r 



7(0, 7 ) -7(1, 7 ) 



7/(1,7 + 1) 
7(0,7) 1 da|k| 



u dr 7/(1, 7 + 1) a|k| dr 



(45) 



The force from the wave pressure is outward on average if the 
average over a pulsation period of expression (45) is larger than 
0. The functions / depend on the wave profile function / and 
are defined by (Koninx & Pijpers, 1992) : 



I(n,j) 



1 

27 



dff- 



(«|k|/0 n 
(l + a|k|/')-> 



(46) 



The integral 1(1, 7+I) is less than zero for realistic profile func- 
tions /, such as a saw tooth profile and a sinusoidal profile, and 
for all values of a Ik I. 



2 (l-l(0, 7 )-l(l, 7 ))/ 7 I(1, 7 +1) 
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Fig. 7. The two terms involving the profile function integrals in equa- 
tion (45) as a function of wave amplitude a|k|. This is for a saw tooth 
profile funtion /. 



The two terms involving the integrals / for the geometrical 
factor 1/r and the velocity gradient term of equation (45) re- 
spectively are shown in figure 7 for a saw tooth profile / and for 
two values of the adiabatic exponent 7. The same factors for a 
sinusoidal profile are not shown because they are virtually iden- 
tical. It is clear that these terms are not negative for any value of 
the wave amplitude and therefore will also not be negative after 
averaging. The force from the sound waves can therefore point 
inward only if the average over a pulsation period of the log- 
arithmic gradient of the wave amplitude is positive and larger 
than the first two terms in equation (45). 

Although the waves damp as they propagate in the direction 
downstream from the shock this does not necessarily mean that 
the gradient of the amplitude a |k| must be positive everywhere. 
If the wind is exactly periodic, as considered in the previous 
section, a shock wave that passes must conserve the sign of this 
gradient. If a region in the wind has a negative gradient of the 
sound wave amplitude this gradient thus remains negative. The 
amplification of the sound waves by the shock is proportional to 
the sound wave action density. The constant of proportionality 
depends only on the amplitude of the pulsational shock. So if 
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to the passage of the shocks. The factor 7a is the ratio of the 
wave action densities behind the shock and in front of the shock 
and it can be calculated from (24) for an adiabatic shock and 
(26) for an isothermal shock. Even for a strictly periodic pulsa- 
tion cycle the strength of the shocks can be a function of radius, 
and therefore so is Ta- The reduced phase angle a(r) is intro- 
duced because the shocks pass two different points r and r' in 
the wind at different phases of the pulsation cycle. Note that 
the occasional input of wave action from waves of type B/C 
below the maximum radius rf to which they can propagate is 
ignored in (34). Only waves of type A are taken into considera- 
tion. The general transport equation for the wave action density 
is more complicated to solve so for clarity the linearized equa- 
tion is used here instead. As was shown by Pijpers and Koninx 
(1992) this means that the results are a reasonable approxima- 
tion as long as the velocity amplitude of the waves is less than 
about half the sound speed. 

The only trivial solution of (34) is A = 0. That solution 
is unstable and the effect of sound waves on a stellar wind tra- 
versed by shocks needs to be modelled explicitly using time 
dependent hydrodynamical simulations with the appropriate 
sound wave pressure term added to the transport equations of 
momentum and of energy, and with (34) as an extra equation to 
solve at every time step. However it is useful to consider sub- 
stituting a solution for the wave action density of the form : 



A = A (t) 



exp 



(r - ro) 



Leg 



(35) 



Here V sw o is the sound wave propagation velocity at radius ro 
from the star, averaged over a period according to (8). This ex- 
pression (35) is then substituted in (34) and all the terms are 
averaged over a stellar pulsation period in the sense of equation 
(8). This leads to the following equation : 



1 ( Ao(t + P) 
p n V Mt) 



^eff 



-1 + 

~L 



r () dL eS 



L e g dr 
In (?a) 



(36) 



+ 



P 



If it is assumed that A is a strictly periodic function : 
A (t + P) = A (t) 



(37) 



then the first term in (36) is identical to 0. It is useful to rewrite 
(36) using the following dimensionless quantities : 



A = L / Rtf 

A e ff = Leg/ i?* 

£ = r/R* 
£o = r /R* 



(38) 



where is the stellar radius. The differential equation (36) 
then becomes : 



1 

A e ff 



■g — go ^ Aeff 



R* In (7a) 

PV SW 



(39) 



Note that this equation does not apply at the sonic point where 
V sv/ = 0. At this point the flow does not transport wave ac- 
tion density and the formulation of the weak damping used 
in equation (34) yields no damping. There is therefore appar- 
ently only an increase in wave action density due to the peri- 
odic impulsive increase due to passing shock waves. This is of 
course un-physical. Strong damping of non-linear sound waves 
is not taken into account here. The non-linear damping terms 
will limit the wave action density at the sonic point to a finite 
value. 

Equation (39) is a form of a Bernouilli differential equation. 
The solution is : 



A e ff 



A 

Co 



1 R* In (?a) 



PV„ 



d£' 



(40) 



The constant of integration is set equal to because A e ff is as- 
sumed to be finite at £ = £o. The term between brackets in (40) 
can easily be recognized as an average of the integrand over the 
region between dimensionless radii £o and £. Note that it is al- 
lowed to have A be a function of radius. The following notations 
for these averages can be introduced : 



1 

1 



i r i 



g-W A 

Co 



d£' 



(41) 



1 [RAn(T A ), cl 
d£ 



£ - go J PV S 

Co 



With these definition the expression for A e ff becomes : 



A 



eff 



(A,) - (A) 



(42) 



If (A*) < (i.e. if V sw < 0) the effective damping length is 
< A e ff < (A). This is the situation in which the sound waves 
of type A can 'drain' back into the star. If < (A*) < (A) the 
effective damping length is < 0. This means that the wave ac- 
tion density locally increases exponentially with radius, rather 
than decreases. For all (A* } > (A) the effective damping length 
A e ff is larger than the 'intrinsic' damping length (A). 

(A*) is not constant as a function of radius, but a rough es- 
timate of the order of magnitude of (A*) can be obtained by 
substituting typical values for cool giants in equation (41) : 



(A* 



0.124 



(T 3 / 100 d)(U sw /lkm/s) 
(7?*/100 R )ln(T A ) 



(43) 



Using again the results from the 'standard model' of Bowen 
(1988) with P = 350 d, 7?* = 270 R Q , and an estimated 
V sw « 0.1 km/s and In 7a = 1, the value of (A*) « 0.02. This 
is roughly equal to the static density scale height of Bowen's 
'standard model'. 

Equating the two different length scales yields the maxi- 
mum effective damping length. As far as the author is aware 
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4 S) 



(1 - Mi) 2 
1 + 2M 2 + 1 /Mf 



+ 1 



1 + 



1 

M 



(27) 



and the ratio of wave action densities by : 



AA ( 7 +l) M 2 + 1 

Ai) B ~ 2 7 Mf - (j - 1) Mi - 1 




(28) 



For transition through an isothermal shock the pressure pertur- 
bation ratio is given by : 



i-\ = M f (l ~ Ml) 
pi 1 1 



[M_ 



(7-D -i 

(2- 7 )J 



(1 + Mi) [Mf + 



Mi 



(7-D 

B Li "^ ' (2-7) T (2-7) 

and the ratio of wave action densities by : 



(29) 



•4i 



1 Mi +1 



Mj 3 Mi 



1 




(30) 



The transmission coefficient as a function of shock strength is 
shown in figure 5. It is clear that in the limit where there is 
no discontinuity and the ratio of the pressure approaches unity 
it should be impossible to convert waves of type B into waves 
of type D. The transmission coefficient indeed decreases to 
in this limit. For very strong shocks the Mach number of the 
material flowing into the shock is very high. It should therefore 
matter little whether a sound wave is of type A or type B since 
the propagation velocity with respect to the shock is almost the 
same. Comparing figures 3 and 5 does indeed show that the 
transmission coefficients for the two wave types are of the same 
order of magnitude. For isothermal shocks there is also a zero 
transmission coefficient for waves of type B for : 



P2 
Pi 



1 



1 



7 



(31) 



It is possible that this is due ignoring the perturbation of the ra- 
diative term in the energy equation at the shock for sound waves 
passing an isothermal shock. If the shock front is radiative and 
the presence of sound waves is properly taken into account in 
the radiation term the ratio of the wave action densities may 
remain finite for all shock strengths larger than unity. 

For waves of type C the transmission coefficient is zero and 
the reflection coefficient can be shown to be (Landau & Lifshitz, 
1987) : 



At _ 1 + M 2 (I- 2M 2 + 1 /Ml 
Ai ~ 1 - M 2 V 1 + 2M 2 + 1 /Mf 



(32) 



This relation holds for reflection at adiabatic shocks and isother- 
mal shocks, but one should keep in mind that the relation be- 
tween M 2 and Mi is different for the two cases. 
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Fig. 6. The ratio of the wave action density of reflected and incident 
waves as a function of shock strength, for case C of fig. 1. The full 
line is for an adiabatic shock, the dashed line for an isothermal shock 

In the limit of a vanishingly small discontinuity the reflec- 
tion coefficient should become zero since there is no disconti- 
nuity to interfere with the sound waves any more. For isother- 
mal shocks the ratio of the densities becomes infinite for infinite 
shock strength. This must have the consequence that waves are 
reflected perfectly and therefore the reflection coefficient must 
approach 1. For adiabatic waves the density ratio remains fi- 
nite and therefore the reflection coefficient remains less than 1. 
In figure 6 the reflection coefficient is shown as a function of 
shock strength. It shows the correct limiting behaviour. 



4. Implications for winds driven by sound wave pressure 

4.1. A periodic solution for A 

From the previous sections it is clear that the presence of strong 
shocks in the atmospheres and winds of red giants does not 
preclude the presence of sound waves. In the modelling of the 
winds of such stars it is therefore important to take both effects 
into account. The transport equation for the wave action density 
of sound waves is (cf Koninx & Pijpers, 1992) : 



dpA 
dt 



+ V • [pAu + B] = pT 



(33) 



Here B is a flux term and T collects all the source terms. For 
linear sound waves this equation reduces to : 



OA 
~dt 



+ V-(tw4) = 
v sw A 



L 



Aln(T A )8 | t 
n = 1, 2, 3, 



afr), „ 
(n + -^)P 

2-7T 



(34) 



Here i; sw is the (instantaneous) propagation speed of the sound 
waves in the rest frame of the star. The second term on the right 
hand side of (34) represents the weak damping of linear sound 
waves in between the strong shocks. The first term on the right 
hand side represent the periodic sound wave amplification due 
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Fig. 3. The ratio of the wave action density in front of and behind the 
shock as a function of shock strength, for case A of fig. 1 . The full line 
is for an adiabatic shock, the dashed line for an isothermal shock 
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Fig. 4. The ratio of the wave energy density behind and in front of 
the shock as a function of shock strength, for case A of fig. 1 . The full 
line is for an adiabatic shock, the dashed line for an isothermal shock 



Here the factor 1/2 is derived assuming a sinusoidal wave pro- 
file. For a saw-tooth profile of the sound waves this factor 
should be 1/3. Using equation (21) the ratio of the wave action 
density in front of the shock and behind the shock can easily be 
calculated : 



( 7 +l) 



l + M 2 



2yM? 



(7 - 1) 1 + Mi 




(24) 



Equation (24) is the wave action density transmission co- 
efficient appropriate for waves of type A of figure 1 passing 
through an adiabatic shock. 

For a strongly radiative shock the perturbation of the radia- 
tive term in the energy equation must also be taken into account. 
If it is assumed that this term is negligible for (nearly) isother- 
mal shocks the ratio of the wave action densities can easily be 
calculated. The pressure perturbation ratio becomes : 



4 S) 



= Mf(l+Mi) 



[Mi + 



(7-D l 
(2-7) 1 



[M? + 



Mi 



(7-D l 
(2- 7 ) ' (2- 7 )J 



(25) 



The wave action density ratio becomes : 



-4i 




(26) 



The reflection coefficient is zero for both adiabatic shocks and 
isothermal shocks. 

In figure 3 the transmission coefficient as a function of 
shock strength is shown for an adiabatic shock and an isother- 
mal shock. The value of 7 is taken 5/3 in both cases. This means 
that away from the shock the gas is assumed always to behave 
as an ideal mono-atomic gas. It is clear that if there is no pres- 
sure discontinuity nothing happens to the sound waves and the 
ratio of the wave action density behind and in front of the shock 
front is unity. Because the density ratio and the velocity ratio are 



higher for isothermal shocks than for adiabatic shocks, isother- 
mal shocks are more efficient at amplifying sound waves. 

The wave action density is the proper quantity to measure 
the effect of shocks on sound waves since the wave action den- 
sity is conserved in the absence of mechanisms generating or 
dissipating waves. However, it is easier to interpret the effect 
of shocks in terms of the energy density in the sound waves. 
This can easily be calculated from the definition of the wave 
action density for linear sound waves (22). The result for the 
waves of type A is shown in figure 4. 
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Fig. 5. The ratio of the wave action density behind and in front of the 
shock as a function of shock strength, for case B of fig. 1 . The full line 
is for an adiabatic shock, the dashed line for an isothermal shock 



For waves of type B the reflection coefficient is also equal 
to 0. The transmission coefficient for type B waves being over- 
taken by the shock differs slightly from the one for type A 
waves. For transition through an adiabatic shock the pressure 
perturbation ratio is given by : 
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fects would put the simulations more firmly into the region of 
fig. 2 where sound waves are dragged out into the wind. 

It is important to note that both in the linear and in the non- 
linear wave pressure tensor (eqs. (1) and (4)) the sign of the 
wave vector cancels out if it is assumed to be purely in the radial 
direction. It is therefore immaterial that the wavevector points 
in the direction of the star in the comoving frame. The only re- 
quirement is that the wave energy density is not equal to and 
this requirement is satisfied. 

Summarizing : up to a point n in the wind waves of all types 
occur and waves of type A/D are generated. Beyond this point 
only waves of type A/D can survive. The strong shock waves 
generated by the pulsation can drag sound waves of type A/D 
superposed on the flow with them so that these will move out 
into the wind. The condition for this in essence is that the aver- 
age over a full pulsation cycle of the outflow velocity is larger 
than the average over a full pulsation cycle of the sound speed. 
The results of numerical simulations to date (e.g. Wood, 1979 ; 
Bowen, 1988) do allow this to occur. It is immaterial for the 
wave pressure that only those sound waves with a wavevector 
pointing towards the star in the fluid frame (type A/D) can sur- 
vive the passage of the strong shock waves. 



3. The amplification of sound waves 

On the passage through a shock the amplitude of the sound 
wave increases. The proof of this is readily available (Landau 
& Lifshitz, 1959) but is repeated here for convenience. The 
starting point is the solution of the equations of conservation 
of mass, of conservation of momentum, and of conservation of 
energy, where the plane sound wave is assumed to impinge on 
the shock along the normal. The coordinate frame is chosen to 
comove with the shock (see fig. 1). These equations are : 



p\V\ = p 2 V 2 

Pi + PlV i = P2 + P2v\ 

1 2 1 2 
Wi H V, = W2 H l>i 

2 1 2 2 



(14) 



If the velocity amplitude of the sound wave impinging on the 
shock is 6v\ the shock wave will start to oscillate with an ampli- 
tude hu and the equations above must be modified. In a linear 
perturbation analysis the equations become 



v\bp\ + p\{bv\ — 8u) = V2&P2 + P2(6v2 — 6u) 
8pi + v\bp\ + 1p\v\(bv\ — 8u) = 

8p2 + v\bp2 + 2p2V2(6v2 — 6u) 

bw\ + v\{bv\ — 8u) = 8w2 + V2(6v2 — 6u) 
The incident sound wave is assumed to be adiabatic : 



(15) 



Just behind the shock there are two types of disturbances. There 
is an entropy wave which is stationary with respect to the fluid. 
There is also a sound wave. The disturbances behind the shock 
are the sum of the two where (S) refers to the sound wave and 
(e) refers to the entropy wave : 

8 P2 = bpf + bpf 
6v 2 = Svf + Sv^ 



02 = 



6w2 = Swf + Sw. 



2 
(e) 

(e) 
2 



(17) 



The relations between the quantities in the sound wave are : 



8 s 



(S) 



8n 



(S) 





CS2 



P2 



£ (S) 



of 



CS2P2 



(18) 



8w 



(S) 



of 



P2 



and for the entropy wave the following holds : 






6pf 

f. (e) 

8vy 



8w 



(e) 



T 2 bsf 



6pf 



(19) 



-S2; 



(7 - l)/°2 

The ratio of the amplitudes of the pressure perturbation of the 
sound waves Spf jbp\ can be determined from the set of linear 
equations (15). It is : 

2M 1 M 2 2 (M 1 2 - 1) - (1 + M0[M 2 + ^Ijj] 
2M 2 (M 2 - 1) - (1 + M 2 )[M? + =^Z] 




Mi + 1 

x 

M 2 + l 

(l + Mi) 2 
\+2M 2 + \/Ml 



1 



7+1 



1 

Mi" 



(20) 



A more general and somewhat more lengthy derivation is given 
in the second edition of 'Fluid Mechanics' of Landau & Lifshitz 
(1987). The velocity perturbation is : 



c (S) 



2 + ( 7 -1)M? 



L2 7 M 2 - ( 7 - 1) 



1/2 



1 



of 



Mi 6pi 



(21) 



For the purposes of sound wave driven winds the quantity of 
interest is the wave action density A which is conserved in the 
absence of wave excitation and dissipation. For small amplitude 
sound waves it is : 



A 



(22) 



8si 
8vi 



8w 





csi 



Spi = 

Pi csipi 

8 pi 
Pi 



(16) 



Here £" w is the wave energy and Q = \w — u • k| is the intrinsic 
wave frequency. For waves propagating in the same direction 
as the flow this reduces to : 



A 



1 . 2^ + c s) 
-pbv 

2 cues 



(23) 
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Bowen (1988) show this behaviour for the velocity also. The av- 
erage velocity in the Eulerian reference frame of a sound wave 
in this period P on emerging from the shock is : 



1 

P .IV 



v(t) - c s (t) dt 



(8) 



If V sv/ < a sound wave starting out in the wind can move 
further towards the star than it gets dragged outwards by the 
mean flow between the strong shocks. In this case the passage of 
one of the strong shocks absorbs the sound waves moving in the 
same direction as the shock and the sound waves propagating 
in the other direction can 'drain' back into the star. If T/ sw > 
a sound wave can only travel part of the way towards the next 
discontinuity and it is dragged away from the star by the mean 
flow. Even if a sound wave starts out in the stellar photosphere 
with a wave vector pointing towards the star the flow between 
the strong shocks will drag it out if V sw > 0. Combining the 
equations above yields : 



V sw = CS1; 



1 



2tt 



M(4>)- 



csi 



d<f> 



(9) 



Here <f> = 2irt / P is the phase of the pulsation cycle. All ve- 
locities are now dimensionless Mach numbers since they are 
divided by the sound speed. The shock speed is divided by the 
sound speed in front of the shock csi. The Mach number of the 
flow M(4>) is : 



M(<f>) = M s 



Mi 
~S~ 



Mi 1 



(10) 



Here S is the ratio of densities behind and in front of the shock. 
Making use of the continuity equation at the discontinuity there 
is a relation between v\ and V2 : 



pi 
Pi 



V2 



(7 + 1 



(7 



l)M? + 2 



S 



(11) 



which holds if the shock waves are adiabatic. If the shocks ra- 
diate, which is more realistic, the law of conservation of mass 
still holds but S is larger. 

In Miras the shocks in the outer atmosphere and wind are 
radiative. The temperature variations are largest for adiabatic 
waves in which case : 



72 = (2jMf - ( 7 - 1)) ((7 - l)Mf + 2) 
Ti ' (y + l) 2 M? 



(12) 



For radiative shocks the temperature ratio is smaller. In the lim- 
iting case of isothermal shocks the ratio is of course unity. For 
isothermal shocks and a flow that is also isothermal the ratio 
csW/csi = 1- For adiabatic shocks and flow the ratio can 
be calculated from the equation of conservation of energy and 
(10): 



'cswr 

, csi , 



1 



(7-1) \Mf - (M(4>) - MsT +1 



(13) 
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Fig. 2. The minimum Mach number of the shock Ms for which sound 
waves are dragged out into the wind. The full line is for a flow with 
adiabatic shocks, the dashed line for a flow with isothermal shocks. 
The dash-dotted line is the minimum Ms for which the average fluid 
velocity is outward for a flow with adiabatic shocks. The dotted line is 
the same for a flow with isothermal shocks 



The condition for V sw leads to a lower limit for Ms so that 
sound waves cannot drain back to the star. This limit can be 
calculated by carrying out the integration in (9) using (10) and 
also the equation of conservation of energy to determine cs(</>). 

The result is shown in figure 2. Also shown in 2 is the 
minimum shock Mach number for which the average fluid 
flow is outward. If Ms is less than the values shown with the 
dash-dotted curve, for adiabatic flow, and the dotted curve, for 
isothermal flow, the average fluid flow is towards the star which 
describes accretion instead of mass loss. There is therefore only 
a small range of shock speeds for which the sound waves can 
'drain' back into the star and this range becomes smaller as the 
shocks radiate more. 

This proves that waves propagating 'towards' the star in the 
frame in which the shock is stationary can actually be dragged 
away from the star by the mean flow. Therefore in general sound 
waves do propagate away from the star albeit 'backwards'. As 
a consequence the wave energy density is not equal to in the 
stellar wind. 

As an example one can use the results of Bowen (1988) to 
estimate the shock strength at around 2 stellar radii and com- 
pare the minimum required shock speed with the actual speed 
of the shock. In Bowen's model without dust (his fig. 5) the 
shock strength in terms of a pressure ratio is around 15. The 
propagation speed of the shock is around 15 km/s which is 
close to the minimum speed for isothermal shocks above which 
sound waves are dragged out by the flow. This makes it very 
hard to conclude whether or not the simulations yield V sv/ > 0. 
In Bowen's figure 5 the shocks are clearly not isothermal but 
the radiative cooling law was necessarily primitive. A more 
complete treatment of the cooling would probably make the 
shocks more radiative and therefore more closely isothermal. 
This would also reduce the pressure ratio at the shock. Both ef- 
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be valid. For an isothermal saw-tooth wave this means that 
the pressure ratio at each discontinuity should be less than 
3/2 + \/2s/5 k. 2.6. For adiabatic saw-tooth waves in a gas 
with a polytropic equation of state the pressure jump should be 
less than 1 + 1 /4j(j + 1) + 1 /4^( 7 + l) 2 + 16. This is approx- 
imately 3.31 if the polytropic index j = 5/3. 

2. The propagation of sound waves 

The propagation of sound waves in an outflow traversed by pe- 
riodic strong shocks can be described most conveniently in two 
steps. The first step concerns the behaviour of the wave at the 
discontinuity. The second concerns the propagation through the 
continuous part of the wind between the shock discontinuities. 

In the coordinate frame where the strong shock wave gener- 
ated by the pulsation is at rest the inflow into the shock is super- 
sonic and the outflow is subsonic as long as (d 2 V/dp 2 ) s > 
where V = 1/p is the specific volume of the gas (cf. Landau & 
Lifshitz, 1959). 




Fig. 1. The velocities and densities in the coordinate frame in which 
the shock is at rest. Note the sign definitions so that the values of vi 
and V2 are positive 

A sound wave with a wave vector pointing away from the 
discontinuity in the upwind direction (case B in fig. 1) will 
therefore be overtaken by the shock. Of course a sound wave 
with a wave vector pointing towards the discontinuity from the 
upwind direction (case A) will reach the shock as well. Both 
types of sound waves will pass through the shock and can then 
propagate away from it with the flow in the downwind direc- 
tion. The amplitude of the sound waves will be increased after 
passage through the shock which is shown in the next section. 
Sound waves can also propagate towards the discontinuity in 
the downwindpart of the flow (case C of fig. 1) since this flow is 
subsonic. Upon reaching the discontinuity these waves must re- 
linquish their energy to the shock wave or be reflected since they 
cannot propagate into a flow that is supersonic in the direction 
opposite to their wave vector. It can be shown that such waves 
are partially reflected (Landau & Lifshitz, 1987). Naively the 



cases B and C are the ones encountered when shocks and sound 
waves both propagate away from the star. 

It is important at this point to realize that the analysis pre- 
sented above is done in the frame comoving with the shock. The 
frame in which the star is at rest is the coordinate frame in which 
the equations for the wind are derived. Therefore it is important 
to know the direction of propagation in this coordinate frame. It 
is clear from the above that after the passage of just one shock 
through a given parcel of gas only the sound waves of type A 
(or equivalently D) survive. This does give waves of type B (or 
equivalently C) some time to travel partway out from the zone 
where the sound is generated. Determination of the maximum 
radius at which waves of type B/C can occur is not trivial. The 
sound waves are most likely to be generated at the top of the 
convectively unstable layer at radius r c . If rs is the distance the 
strong shock wave travels in a time t after passing r c then it 
satisfies : 



drs 
rs = 



vs(rs,t) 
~ at t = 



(5) 



Here vs is the velocity of the shock wave in the rest frame of 
the star. For the sound waves of type B/C emitted from the top 
of the convective layer at some phase delay a after the passage 
of a shock, the distance a wave front can travel is : 



dr 



dt 



= v(r sw ,t) + c s (r sw ,/) 
a 



(6) 



at/ 



-P 



iir 



Here v is the velocity of the fluid flow in the rest frame of the 
star, cs the sound speed, and P the period between passages of 
the strong shock waves. The sound waves are overtaken by the 
next shock when r sw = rs = r„. The maximum distance that 
sound waves can travel out from the star is then n = max r a 
as a function of a. Although it is possible in principle that n 
is infinite, it is likely that beyond some finite radius all the gas 
will have passed at least once through a shock front. The region 
where waves of type A/D are generated, either directly from 
turbulent convection or through shock processing of waves of 
type B/C, extends up to this radius rs . Beyond that radius only 
waves of type A/D survive. In the rest of this section therefore 
only the case A of figure 1 is considered. 

In the frame in which the star is at rest the velocity of the 
flow emerging from the discontinuity at time / = at a fixed 
position with respect to the star is : 



v(t) = Vs — V2 + (l>2 — v l) 



P 



(7) 



The signs have been chosen such that all velocities on the right 
hand side in equation (7) have positive values. On the left hand 
side v(t) is positive for material moving away from the star. It 
is assumed that the velocity is linear as a function of distance 
behind the shock and therefore also as a function of time. This 
linearity is a very good approximation for strong shocks (cf. 
Landau & Lifshitz, 1959). The hydrodynamical simulations of 
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Abstract. It has been claimed that in stellar winds traversed by 
strong shocks the mechanism for driving the wind by sound 
wave pressure cannot operate because sound waves cannot 
propagate past the shocks. It is shown here that sound waves 
can propagate through shocks in one direction and that this is a 
sufficient condition for the sound wave pressure mechanism to 
work. A strong shock amplifies a sound wave passing through 
it and can drag the sound wave away from the star. It is imma- 
terial for the sound wave pressure gradient that the sound wave 
vector points towards the star. Since the strong shocks drag the 
sound waves away, the star itself is the source for the sound 
waves propagating towards it. 

Key words: Shock waves - Stars : oscillations - Stars : mass 
loss 



1. Introduction 

The effect of the presence of sound waves on the winds from 
stars can be modelled by introducing a wave pressure (Pijpers 
& Hearn, 1989). The gradient of this pressure acts in the same 
way that the gradient of the gas pressure does. For a linear sound 
wave which has a small amplitude the pressure can be calcu- 
lated explicitly without a knowledge of the background veloc- 
ity as long as this is smooth. The appropriate expressions are 
derived by Pijpers and Hearn (1989) on the basis of the linear 
perturbation formulations for the propagation of sound waves 
in continuous moving media of Bretherton (1970). 

The linear theory of sound wave driven winds has been used 
to model the wind of cool giants (Pijpers and Habing, 1989). 
It has also been used to model to outflow from the M super- 
giant VX Sgr (Pijpers, 1990). Objections have been raised (e.g. 
Wood, 1990) to the application of this model to the Mira class 
of variable M giants. The reasoning is that the stellar pulsation 
of Mira variables generates strong shocks which propagate out 
into the wind. These strong shocks would overtake and absorb 
the sound waves so that the pressure from sound waves would 
be effectively cancelled. The purpose of this paper is to demon- 
strate that this is incorrect. 



The expression for the wave pressure tensor in the linear 
regime is : 

*U- = -r^ 2 | k | 2 (^ + ^T-D%) (1) 

Here k is the wave vector of the sound wave and Sij is the Kro- 
necker delta, cs is the speed of sound defined by 



where p is the ordinary gas pressure and p the gas density, p is 
the Lagrangian mean density of the flow, a |k| is the amplitude of 
the displacement in the sound wave, j is the adiabatic index. If 
the waves become non-linear the approximations of this linear 
perturbation method are no longer valid. 

It is still possible to derive a pressure tensor for non linear 
waves as long as the velocity amplitude of the sound wave is 
less than the sound speed. To do this one can follow the gen- 
eralized Lagrangian mean formalism (GLM) presented by An- 
drews & Maclntyre (1978a, 1978b). Koninx & Pijpers (1992) 
show that if the disturbance profile of the sound waves is written 
as : 

( = a±-f(k-x-iot-6) (3) 
l k l 

the general wave pressure tensor takes the form : 

i r 1 ' ( % + fe - %H«MT\ 

<-r r l a+ wf)' j M <4) 

/' is the derivative of the profile function /. The limitation to 
the velocity amplitude of the sound wave that is less than the 
sound speed means that a|k| < 1. 

Note that the profile function / gives considerable free- 
dom of the shape of the sound wave profile. The waves are not 
required to be sinusoidal. It is well known that sound waves 
in stellar atmospheres very quickly deform into a saw-tooth 
shape with small discontinuities. As long as the velocity am- 
plitude of the saw-tooth is less than the sound speed the con- 
ditions are met for the non-linear wave pressure formalism to 



